Quantum instability in a dc-SQUID with strongly asymmetric dynamical parameters 
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A classical system cannot escape out of a metastable state at zero temperature. However, a 
composite system made from both classical and quantum degrees of freedom may drag itself out 
of the metastable state by a sequential process. The sequence starts with the tunneling of the 
quantum component which then triggers a distortion of the trapping potential holding the classical 
part. Provided this distortion is large enough to turn the metastable state into an unstable one, the 
classical component can escape. This process reminds of the famous baron Miinchhausen who told 
the story of rescuing himself from sinking in a swamp by pulling himself up by his own hair — we 
thus term this decay the 'Miinchhausen effect'. We show that such a composite system can be 
conveniently studied and implemented in a dc-SQUID featuring asymmetric dynamical parameters. 
We determine the dynamical phase diagram of this system for various choices of junction parameters 
and system preparations. 

PACS numbers: 85.25.Dq, 74.50.+r 



I. INTRODUCTION 

Consider a classical object (a degree of freedom) 
trapped in a metastable potential minimum; no decay out 
of this metastable state is possible at low temperatures, 
where thermal activation over the barrier is exponentially 
suppressed. However, if the classical object is a compos- 
ite one, with a quantum degree of freedom coupled to 
the classical one, then the quantum object may tunnel 
out of the metastable minimum and exert a pulling force 
on the classical object. Once the latter is large enough 
to completely suppress the trapping barrier, the classical 
object is able to leave the potential well — hence a clas- 
sical object may escape from a metastable state even at 
zero temperature if helped by a coupled quantum degree 
of freedom. 

The above situation can be realized in a dc-SQUID 
(Superconducting Quantum Interference Device) featur- 
ing asymmetric dynamical parameters; i.e., with two 
Josephson junctions of equal critical currents J c but 
strongly different (shunt) capacitances C and (shunt) re- 
sistances R, see Fig. [T] choosing large and small parame- 
ters C and 1/R for the two junctions allows to place one 
of the junctions in the 'classical' and the other into the 
quantum domain. The tunneling of the quantum degree 
of freedom entails a distortion of the trapping potential of 
the classical junction, which might be sufficiently large to 
transform the metastable state of the classical junction 
into an unstable one. The appearance of this complex 
decay channel depends critically on the applied bias cur- 
rent J and the SQUID's loop inductance L coupling the 
two junctions. 

The gauge invariant phase differences ifi, i = 1,2, 
across the two Josephson junctions^ define our dynami- 
cal degrees of freedom: Assuming equal critical currents 
J c , the junctions' potential energies V; = Ej[l — cos ifi], 
i = 1,2, involve the Josephson energy Ej = <& J c /2-kc 
(with $o = hc/2e the flux unit, e and c denote the 
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FIG. 1: Schematics of the dynamically asymmetric dc- 
SQUID. Two Josephson junctions are integrated in a current 
(J) biased superconducting loop with inductance L. The two 
junctions feature equal critical currents J c but strongly asym- 
metric (shunt) capacitances d and (shunt) resistances Ri. 
We assume that d and Ri are chosen such that quantum ef- 
fects are present for junction 2 but are strongly suppressed 
for junction 1, hence Ci 3> C2 and /or Ri <C i?2- 



unit charge and light velocity)^ Their kinetic energies 
% = (h/2e) 2 Ciifl [2 are determined by the junction ca- 
pacitances Ci playing the role of effective masses (the 
relevant energy scale is given by the charging energy 
E c .i = e 2 /2Ci) — a dynamically asymmetric SQUID with 
one large and one small junction capacitance then pro- 
vides us with the desired classical and quantum degrees 
of freedom (we choose C\ 3> C2; additional normal re- 
sistances Ri introduce a dissipative dynamical compo- 
nent, see below). The coupling of the two junctions 
via the loop inductance L produces the interaction en- 
ergy Vint = [($ /27r) 2 /2L](</3 2 - </?i) 2 involving the rela- 
tive coordinate (^2 — fi), whereas the external driving 
current J couples to the absolute coordinate, Vdrivo = 
Ej(J/2J c )(ipi + (p2)- While large-capacitance ('classical') 
junctions are easily fabricated, small ('quantum') junc- 
tions are more difficult to realize. Nevertheless, experi- 
mental techniques to fabricate small junctions are avail- 
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able today and their quantum behavior in the form of 
quantum tunneling^ 3 -^^ quantized energy levels^ and 
even quantum coherenc e) 9 ' 10 ' 11 ' 12 ' 13 ' 14 ' 15 has been demon- 
strated. 

The quantum decay of the biased, dynamically 
symmetric dc-SQUID has been discussed before both 
experimentally 16 .^ and theoretically? 1 ^ also in the con- 
text of instanton splitting! 19 ' 20 In the SQUID discussed 
here^i the large asymmetry of the dynamical parame- 
ters blocks the tunneling of the ipi degree of freedom; the 
'heavy' junction then remains frozen during the quantum 
decay of the 'light' degree of freedom <p2. A trajectory 
of this kind describes the entry of flux into the SQUID 
loop which rearranges the current flow in the two arms 
in a way as to redirect more current through the heavy 
junction. This increase in current produces an enhanced 
tilt — Jcffi^i in the potential of the heavy junction, which 
then may decay through a classical trajectory, see Fig. [5] 
We call this nontrivial decay sequence the 'Miinchhausen 
decay'. It is the aim of this work to determine the effec- 
tive critical current J C (L _1 ) for which the Miinchhausen 
decay becomes possible, see Figs. 151 and 

In the following, we define our system in full detail, in- 
cluding a dissipative component in the junction dynamics 
(Sec. |ll]). Sections UTTl and [TVl are devoted to the deriva- 
tion of the effective critical currents for the various cases 
with junctions governed by massive or dissipative dynam- 
ics. In Section [V] we add remarks concerning the experi- 
mental realization of the system described here. Finally, 
we draw conclusions in Section IVIl 



II. SETUP AND MODEL 

Within the resistively and capacitively shunted junc- 
tion (RCSJ) model (at T — 0), the classical dynamics of 
the two phase differences tpi and ip 2 is governed by the 
equations of motion 

u >o ! i<Pi + m l j>% = -d<piV{<P\,<P2), (1) 

with the plasma frequency h 2 uj^ i = &EjE c ^ of an un- 
biased single junction and the damping coefficients r\i — 
^o/2TrcJ c Ri (Ri denote the normal ohmic junction resis- 
tances). The potential (see Fig. [5]) is given by 

v(<pi, 1P2) = 1 - cost^i + 1 - COS</3 2 

k 

- j((p l +(p 2 ) + -((fix- ip 2 f, (2) 

with the dimensionless current j = J/2J C and the cou- 
pling constant k — $oc/2ttJ c L = 1//3l {@l denotes the 
usual screening parameter of the SQUID). Eqns. fl} and 
<j2j) describe a dc-SQUID with symmetric inductance L in 
a vanishing external magnetic field and driven by a bias 
current J or, equivalently, the massive (mass cx w^l) 
and/or dissipative (r^) dynamics of two harmonically (k) 
coupled particles in a tilted (J) and corrugated (cos ifi) 



potential. Quantum effects of the light junction 2 are 
accounted for via the relevant tunneling and decay pro- 
cesses, see below. 
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FIG. 2: Surface (a) and contour plot (b) of the SQUID poten- 
tial, Eq. ([2]), at bias current j = 0.5 and coupling k = 0.04. 
For an undercritical current j < 1, the symmetric minima 
(<pi = 1P2) are all metastable. The stability of the side-minima 
(<pi 7^ ifi2) depends crucially on the parameters j and k. The 
line in (b) shows the decay path in an overdamped setup: The 
system starts out in a relaxed local ground state near the ini- 
tial minimum at tpi = ip 2 = arcsinj. The 'light' degree of 
freedom if2 then tunnels (dotted line) and the system relaxes 
to the bottom of the next minimum near ip 2 ~ 2n, with a clas- 
sically stable but quantum mechanically metastable ground 
state. Through an additional quantum phase slip, the light 
phase reaches the minimum near <p2 ~ 47r, which is not clas- 
sically stable and henceforth the system can decay along a 
classically allowed path to ipi ~ 2n. The system then has 
turned unstable and enters a resistive state through iteration 
of the last two steps. 

The potential © gives rise to two types of relevant 
frequencies. One is the plasma frequency u> Pt i, the small- 
amplitude frequency in the direction of <pi around a local 
minimum of the potential v((pi,<f2)- With the effective 
potential 

Ves[<pi]((pi) = v(tpi = const., ifi), i ^ I, (3) 

cf. Fig. [3l Up i — u^d 2 .i> e ff(^i)j evaluated at a local 
minimum (pf nn , and depends on the parameters j and 
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k as well as on ff 1 ™. For the heavy junction (junction 
1), can become arbitrarily small upon approaching 
criticality, while for the quantum junction (junction 2) 
u p ,2 becomes small only for j — ► 1 and k — ► 0. The 
other frequency is given by the LC constant of the 'su- 
perwell' in v e s{f2), cf. Fig. [3J and is relevant only in the 
regime and for the quantum junction (junction 2), 

W LC,2= W 0,2^C 2 /LC 2 . 
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FIG. 3: Effective potential v e ff(¥>2) = w(vi = const., 9J2) 
(solid line) and the parabola remaining after dropping cos fi 
(dashed line) for j — 0.5, k — 0.02 and 931 = arcsin j. The bul- 
lets and arrows illustrate the sequential decay of the quantum 
phase if 2 to the ground state for the case of strong damping. 

These characteristic frequencies delineate the regimes 
where the two junctions behave classically and quan- 
tum mechanically, respectively. The relevant physical 
parameters are given by the ratios Ej/E c ^. For the 
heavy junction, we have to make sure that the num- 
ber Ni w E b /huj pA ~ (Ej/E c>1 )(1 - Ji/J c ) 5/4 of states 
in the local well remains large. Choosing a large ratio 
Ej/E Ct i 3> 1 then guarantees that quantum effects (i.e., 
Ni of order unity) become relevant only very close to 
criticality and we choose to ignore them in our discus- 
sion below. 

The system develops the most interesting behavior if 
the quantum degree of freedom develops well localized 
states within all relevant local minima in u e ff(</?2), Fig. 
02 and hence parameters should be chosen such that the 
quasi-classical description applies. On the other hand, 
we require tunneling and coherence effects to manifest 
themselves on reasonably short (measurable) timescales. 
We then choose a ratio Ej/fkvo^ ~ 2, such that the local 
wells in UeffCva) contain a few quasi-classical states each. 

The strength of dissipation can be quantified by the 
dimensionless damping parameters, 



CHLC,2 = ■ 



(4) 
(5) 



Below, we are interested in the two limiting cases of 
strong and weak damping. For a strongly damped quan- 
tum junction with a Pt 2 > 1, the quantum decay of if 2 out 
of a metastable well of ^ e ff (^2) is incoherent^*^ and its 
subsequent relaxation is fast (as compared to the dynam- 
ics of ipi). For weak damping, ar,c,2 -C 1, a p ,i <C 1, the 
kinetic energy stored in the motion of the heavy junction 



has to be accounted for; in addition, the finite lifetime of 
the quantum states of the light junction due to the resid- 
ual dissipation has to be considered, see the discussion in 

sec. inn 



III. STRONG DAMPING 

We start by analyzing the situation for strong damp- 
ing, Cf p 2 > 1 and a Pi i 3> 1, see Ref. [25| for an earlier 
short report. We bear in mind an experiment with a dc- 
SQUID characterized by a fixed inductance L oc fc _1 and 
biased with a current j < 1. The task is to determine 
whether the Munchhausen decay can take place; in the 
experiment, the latter manifests itself through the transi- 
tion to a finite voltage state. In the strong damping case, 
no kinetic energy is stored in the system. Furthermore, 
the evolution is not sensitive to the way the current is 
ramped. After current ramping, the system starts out in 
a relaxed state where the phases (fi are localized in the 
diagonal metastable minimum at tpi = if 2 = arcsin j (up 
to an arbitrary multiple of 2tt). 

For sufficiently large j, the quantum degree of freedom 
if 2 undergoes tunneling to a new local minimum nearby 
27rn, while the classical degree of freedom (pi — arcsin j 
remains localized, thus allowing a flux ~ n$o 5 n £ N 
to enter the SQUID loop, cf. Fig. [3] If the resulting 
force on the classical phase f\ is sufficiently large, the 
Munchhausen decay is enabled with a classical decay of 
fi and successive iteration of quantum decay (directed 
along f2, flux entry) and classical relaxation (directed 
mainly along (pi, flux exit), cf. Fig. 03 

The strong coupling at large k keeps the local mini- 
mum at f 2 = arcsin j lower in energy than the adjacent 
local well at (p2 ~ 2-7T for all bias currents j, hence tunnel- 
ing is inhibited and no Munchhausen decay takes place. 
At lower k, a bias current j < 1 can sufficiently lower 
the adjacent well such as to bring both minima to equal 
height. This condition is reached once the minimum of 
the parabola in u e if (V2) (at (f% = arcsin j+j/k) is aligned 
with the midpoint between the two corresponding min- 
ima of — cos</?2 (at (f2 = 7r); i.e., if arcsin j + j/k = n. 
Thus, for k < k+^j) with^S 



J 



arcsin j 



(6) 



the minimum at f>2 ~ 27r is lower and a quantum decay is 
enabled (cf. Fig.[D(b), point C in the diagram). The jump 
of ip 2 by roughly 2-7T then pulls the heavy junction out of 
its minimum and the Munchhausen decay is initiated, cf. 
Figs.[UJa) and (b) associated with the sequence of points 
A, B, and C in Fig.|U[c). The smaller k becomes, i.e., the 
weaker if 2 is bound to ifi, the less current is needed to 
lower the adjacent minimum such as to allow for a decay 
of f 2, hence the positive slope of kj^j). 

Decreasing k too far, however, the pulling force ex- 
erted by the quantum junction may not be sufficient to 
drag the heavy junction out of its minimum. Hence, we 
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FIG. 4: Illustration of the assembly of j c (k) from the segments 
j+i(k), Eq. ©, and j~x(k), Eq. JTTJ. In (c), the relevant 
region of the j-fe-plane is shown. At A, the quantum junction 
(see (b)) has relaxed to the initial minimum (n = 0) and 
the classical junction is localized in a stable minimum (see 
(a)). Increasing the bias current j and approaching B, the 
potential barrier confining the classical phase in v e g((pi) is 
only slightly reduced. As j > jj" 1 (fe) (C), the increase in 
j has deformed the effective potential v e g((p2) as to allow a 
phase slip of (p2 (see (b)). The additional force it exerts on 
the classical junction immediately removes the barrier (see 
(a)) and delocalizes (pi at j c (k) = j^i(k). On the other hand, 
an increase in bias j > (k) from A' to B' at lower k leads 
to a phase slip of tp2 (see (e)), without delocalizing (pi (see 
(d)). A further increase in j triggers a classical decay (see 
(d)) when crossing the critical line at j c (k) = jji(fc) (C). 



have to investigate the shape of the potential v((fx,(f 2 ) 
after the phase slip in (p 2 , i.e., nearby the point (p% > 
arcsinj, (p 2 ~ 2tt, and check whether the barrier against 
a classical decay (mainly along (fx) has disappeared; this 
is identical to the calculation of the critical current of a 
SQUID with a trapped ffux,2L22 

We first determine the position of the minimum by 
solving the equations 

d Vl v((fi, (p 2 ) = sin(^i) - j + k{<px - (p 2 ) = 0, (7) 
d V2 v{(pi,ip 2 ) = sin(<^ 2 ) - j - k((f! - (p 2 ) = 0, (8) 

for ipi > arcsinj and (p 2 ~ 27r. At the critical cou- 
pling k~i(J) the minimum should merge with a saddle 
and define an inflection point along some direction in 
the {(pi, ^2)-plane. The resulting system of equations re- 
quires numerical solution and the result is shown in the 
inset of Fig. [5l However, within the interesting region 
at small coupling k we can find an approximate analyt- 
ical solution: For i<C 1, the side minimum of v(tp\, (p 2 ) 
becomes unstable predominantly along the </?i-direction. 



Thus, the minimum disappears if 

d$ )1 v(<pi,<p2)=co8<pi+k = Q. (9) 

We choose the solution (p\ m tt/2 (we set k = and 
assume < (p\ < 2tt; the other solution (p\ w 3n/2 
cannot solve Eq. ([7]) and is discarded). Inserting (px into 
the sum of Eqs. (J7J) and ([5]) yields the relation 

sin^ 2 «2j-l, (10) 

from which we find <p 2 rs 2n + arcsin[2j — 1] (the other 
solution ip 2 w 3tt — arcsin[2j — 1] is excluded since it 
describes a maximum along the 922-direction). Inserting 
(fx and tp 2 into Eq. (0, we find the condition 



C ' 1U; (3/2)tt + arcsin(2j - 1)' v ; 

For k > k~x(j) no barrier blocks the motion of ipx after 
the phase slip in (p 2 and a classical decay of (fx is enabled, 
thus completing the Munchhausen decay. This scenario 
is described by the sequence A, B, C in Fig.0Jc). For fc < 
k~x(j) force after the phase slips is too small and an 
additional increase in j is necessary to drive the system 
overcritical, as illustrated by the sequence A', B', C in 
Fig. [D^c) and Figs. H{d) and (e). The increase in critical 
current with decreasing coupling defines a negative slope 
for k~xU)- 

We define the effective critical current j c (k) as the 
phase boundary between the stable region (j < j c (k)), 
where the Munchhausen decay is prohibited, and the de- 
localized phase (j > j c {k))- This critical line is assem- 
bled from the segments j^fc), the inverse functions of 
k c -^j), Eqs. © and (fTT|) . respectively: Ramping up the 
current j at large values of fc, we eventually cross j^x(k)- 
The phase slip of <p 2 immediately enables the classical de- 
cay of (fx and the system enters a running state, hence, 
j c (k) — 3cx{k). At lower k, increasing j beyond x (fc) 
triggers a phase slip of ip 2 , but the resulting force is too 
weak to delocalize ip\. A further increase in j is neces- 
sary until, at j = j~x(k), the minimum disappears and 

the Munchhausen decay proceeds, hence, j c (k) = jcx(h)- 
Following the above discussion, the system is stable 
for j < j~x(k). Upon further decreasing k, however, 
more and more side minima in v e s((p 2 ) become acces- 
sible to the quantum junction. From Fig. [3] we notice, 
that the different local minima of v c s((p 2 ) are located 
near ip 2 ~ 27rn. They describe the state where a flux 
~ n$o has entered the loop and we label them by the 
index n. Eq. ^ is then straightforwardly generalized 
to the critical line k^ n {j) describing the entry of the pi- 
th fiuxon. The n-th phase slip of (p 2 occurs when the 
global minimum in v c g((p 2 ) shifts from ip 2 ps 27r(n — 1) to 
ip 2 k 2nn. We can proceed analogously to the derivation 
of Eq. ([6]): Let (/?™ m,rl-1 be the solution of Eqs. ([7]) and 
near (fx ^ arcsinj, <p 2 rj 27r(n — 1). In order to find 
the crossing in the height of the two minima we account 
for the relaxation of (fx and align the minimum of the 
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shifted parabola (if2 = (p I f lm,n ~ + j/k) with the mid- 
point between the two minima of the — cos (/^-potential 
((f2 = (2n — 1)tt), hence ip™ ln ' n ~ + j/k = (2n — l)ir and 



fn 1 \ min.n— 1 

(2n — 1)7T — (f 1 



(12) 



A convenient approximation is made by ignoring the re- 
laxation of ipi, resulting in the expression 



J 



(2n — l)n — arcsinj 



(13) 



Similarly, we can generalize Eq. (|lip to a critical cou- 
pling k~ n (j), determining whether the resulting force at 
given fluxon index n is sufficient to delocalize tpi. We 
then solve Eqs. and © near (pi rj arcsinj, ip2 ~ 27rn, 
proceed as in the derivation of Eq. (fTT|) . and find (for 
small k < 1) 



1-3 



(2n - 1/2)tt + arcsin(2j - 1) 



(14) 
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FIG. 5: Phase diagram of the dynamically asymmetric dc- 
SQUID as a function of bias current j = J/2J C and inductive 
coupling k = &oc/2nLJ c . Here, we assume strong damping 
ot p ,i 2> 1 and a p ,2 > 1. The effective critical current j c (k) 
(solid line) marks the boundary between a localized classical 
junction (lower bias j < j c ) and a delocalized classical junc- 
tion (j > j c ), corresponding to a finite voltage state of the 
SQUID. Branches with negative slope are determined by a 
classical instability (mainly along tp{), while those with pos- 
itive slope are determined by a quantum instability of the 
light junction. For j < j c {k), the dotted lines j^„(fc) mark 
the entry of flux through the quantum junction (the integer 
n approximately quantifies the flux through the ring in the 
stable state); these lines can be measured via monitoring of 
the flux threading the loop. For j > j c (k), the dashed lines 
j~ m (fe) mark the minimum number m of flux units necessary 
to delocalize the classical junction. The inset shows a com- 
parison between the approximate result, Eqs. (fl3|) and (fl4)) 
(solid line) and the exact numerical result (dashed). 

The critical current line j c (k) is constructed from in- 
terchanging segments of j^ n (k) and j c T n (fc), resulting in 
the dynamical phase diagram, Fig. [5] Note that the dif- 
ferent nature of the decay, classical or quantum, associ- 
ated with the two types of critical lines may allow for 



an experimental distinction: Ramping the current past a 
+ -type segment of j c {k) triggers a quantum decay with 
a broad histogram describing multiple measurements. A 
"-type segment of j c (fc) triggers a classical decay with a 
sharp histogram (the quantum decay of ipi needs to have 
occurred already, which limits the ramping speed before 
reaching the critical line). Note that the lines j c f „(fc) are 
detectable throughout all the stable portion of the phase 
diagram, e.g., via a measurement of the flux threading 
the loop (the flux increases by approximately one flux 
unit upon crossing the dotted lines in Fig. [5]) . 

The phase diagram in Fig. [5] shows that the critical 
line j c (fc) approaches the value 1/2 for k — > 0. This is 
understood from analyzing Eqs. (fl"3]) and (fl"4"|) for large 
n and small k, providing the relations kj n (j) ks j/2nn 
and k~ n {j) « (1 — j)/2nn, respectively. Equating the 
two conditions gives j c w 0.5. From a physical point of 
view, this result can be easily explained: The decay of ip2 
proceeds towards the bottom of the parabola in v e s (^2); 
for k — ► 0, the current through the quantum junction 
then approaches zero, J2 oc sin (^2 ~ 0. Consequently, 
all current is redirected through the classical junction 
1, with the effective bias now increased to J\ = 2jJ c . 
Junction 1 thus turns dissipative at j = 1/2, the critical 
current of a single junction but half the critical current 
of the dc-SQUID. 




FIG. 6: Sketch of a time-trace of the flux threading the 
SQUID loop during the decay process of ipi and ip2 at j > 
j c (fe), cf. Fig. [2jb). The flux increases by approximately one 
flux unit each time the quantum junction decays; this pro- 
cess then involves the tunneling rate F2 of ip2, Eq. (|15[) . The 
additional flux leaves the loop again within the time interval 
~ 7/1 as the classical junction relaxes dissipatively to the next 
minimum. 

The dynamical phase diagram, Fig. [SJ can be tested 
experimentally through a measurement of the dc-voltage 
drop across the device. The second Josephson relation 1 
tells that the time-averaged voltage (V)t oc (<pi)t = 
(if2)t- Inspecting Fig.[^b), we understand that the con- 
tinuous iteration of quantum and classical decays of the 
light and heavy phase variables generate the nonzero av- 
erages ((fi)t = (v?2)t for j > jc(k). For the overdamped 
setup as described above, the voltage drop is small as 
strong damping reduces the tunneling rate determining 
(ifi2)t- This is particularly true close to the critical line, 
where the variable (p2 needs to tunnel in the 'flat' part of 
w e ff(^2)- The decay process is faster for weak damping 



6 



and thus generates a larger voltage signal — we discuss 
this situation in the next Section. 

Another characteristic signal of the Munchhausen de- 
cay is the time-trace of the magnetic flux threading the 
SQUID loop during the alternating decay of the two 
phases (cf. Fig.[2jb)), see Fig.[6]for an illustration. Such 
time-traces exhibit two characteristic time scales, one due 
to the tunneling of the quantum junction involving the 
rate T 2 separating subsequent peaks of flux entry into 
the loop; the other is due to the classical relaxation of 
the heavy junction and involves the dissipative time ~ 771 
describing the flux exit and return to the low-flux state 
of the loop. At small k and small effective bias between 
two minima, 

r 2 W) _ ^Q,2^c,2 / _^Q_\ 7/2 (■ _ ■+ ^R Q /R 2 -l Mrs 

with R Q = h/4e 2 M 

IV. WEAK DAMPING 

The behavior of the system for weak damping ollc.i *C 
1, a p ,i <C 1 is quite similar to the one encountered before, 
while the (small) differences strongly depend on the many 
system parameters and their associated timescales. To 
begin with we neglect dissipative effects. Bearing in mind 
the limit of large C1/C2, we can exploit the adiabatic 
separation of fast and slow degrees of freedom; i.e., we 
first consider the fast problem for the quantum junction, 

H 2 ^i{if 2 ) = £;*i(<^ 2 ), (16) 
with the reduced Hamiltonian (cf. Eq. {I)) 

H 2 = -4E C:2 d 2 V2 + £jUeffM(¥>2) (17) 

defined at fixed f\. Eq. p6[) establishes the dynamics 
of the 'light' variable <f 2 and we can find the energies 
Ei(ipi). These act as effective potentials for the 'heavy' 
degree of freedom ipi; e.g., assuming the quantum junc- 
tion to reside in a state |Z), {f 2 \l) = ^>i(f 2 ), the effective 
potential for ipi is given by £/(yi). Typical examples for 
the effective potentials £o(fi) (ground state) and £i(<^i) 
(first excited state) are shown in Fig. [7] 

In the dissipative situation, we could derive the phase 
diagram from simply analyzing the potential v(if\,f 2 )] 
here, instead, we first determine which states I of the 
light junction are relevant and then analyze the dynamics 
of if 1 in the resulting effective potential (e.g., Fig. [T^a)) 
defined via £;(</?i), thus determining whether the heavy 
junction stays localized or enters a running state. Both 
tasks strongly depend on the various timescales involved. 
We start with the timescale of the 'external control', the 
rate at which the bias j is ramped. Whereas the case of 
strong damping excluded the presence of kinetic energy 
in the system, this is no longer the case here. The im- 
portance of the ramping rate then comes about through 



the amount of energy which is transferred to the two de- 
grees of freedom. The change of j translates to a change 
of the potential v{tf\, if 2 ) with a rate of the order of dtj. 
The ramping can be adiabatic, in which case no energy is 
given to the system, or instantaneous, where the amount 
of transferred energy is maximal. Intermediate types of 
ramping will not be discussed. 




tt/2 tt/2 




tt/2 tt/2 
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FIG. 7: Illustration of the effective potential for the 'heavy' 
phase <fi, generated by the two lowest levels of the 'light' 
phase. In (a), the situation for k = 0.15, j = 0.36 is shown. 
The avoided level crossing 2A2 in the spectrum of the light 
junction is reflected in the effective potentials £o{fi) (solid) 
and Ei(fi) (dotted) for the heavy junction. The inset illus- 
trates the 2-dimensional path [tpi, ^ in '°^ 1 (fi)] in v((pi,f2) 
corresponding to eo(fi) (solid, A/B) and ei(fi) (dotted, 
C/D); the upward/downward jump along if 2 to the next mini- 
mum corresponds to the avoided crossing in e(fi). (b) Series 
of £0(^1) (solid) and ei(fi) (dotted) for various values of j 
at k — 0.1. Increasing j shifts the position of the avoided 
crossing in e(yi) to the left. Its transition past <^i" n (be- 
tween j — 0.26 and j = 0.3) corresponds to crossing the line 
ic~i W an( l implies the decay of the quantum junction to the 
next local minimum. Note that in the present case the newly 
formed minimum is stable (a point of type B' in Fig. [4| . 



The adiabatic separation of time scales requires a large 
capacitance ratio C± 3> C 2 . For Lu Pt i <C u> Pr 2 the in- 
travalley motion of the two junctions are separated. The 
(tunneling) motion of the light junction across valleys 
produces a stronger condition: Consider the tunneling 
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splitting A2 in the spectrum of H2 at an avoided level 
crossing, cf. Fig. [TJa), 

A 2 ~ hu Pj2 exp(--,^Ej/4E Ci2 ), (18) 
with 7 depending on the actual shape of the potential, 



1=1 dtp 2 y/v c s (f2 ) ~ E 1 Ej (19) 

J a 

in the quasi-classical approximation. Here, Eq is the 
ground state energy in the well and a and b are the left 
and right classical turning points, respectively, v c g(ip2 — 
a) = v e s(b) = Eq/Ej. The tunneling gap A2 then is rel- 
evant when analyzing the motion of the classical junction 
(ifi) in the effective potential given by £i(if\). If the mo- 
tion of if 1 is sufficiently slow, its passage along an avoided 
crossing in the spectrum of H.2 is adiabatic and ip± follows 
£i(<fi) through the anticrossing (cf. the trajectory A — > B 
in Fig. 0Ja)); i.e., the phase if 2 tunnels to the next ad- 
jacent minimum. On the other hand, if the motion of 
ip 1 is fast, the state of the quantum junction undergoes 
Landau-Zener tunneling (cf. the trajectory A — > D in Fig. 
[TJ^a) ; the light junction then remains trapped in its local 
well) and tpi follows the effective potential £i±i(<fi) after 
the avoided crossing. The probability j»lz for Landau- 
Zcner tunneling at an anticrossing of the two lowest levels 
of H2 is given by 

PLZ = ex H"^MJ' (20) 

where £io(i) = £i(t)~£o(t) and dt is the total time deriva- 
tive. To estimate the rate of change in energy dt£io(t), 
we start from 

~. T7i r / min,n+l\ / min,n\i 

£l w Ej[v(lf U lf 2 )-v(ip 1 ,if 2 )J (21) 

and find the time derivative 

d t Ew(t) = d Vl £iQ(fx 

+Ej [d V2 v((f 1 ,if 2 )\ tp ^m,„+id lpi (f 2 nm ' n+1 

r\ t \ 1 o min.nl ■ 

-d V2 V(<fl,lfi2)\ v ™n,nd Vl (f 2 \ipt 

7 771 / min.n min,n+l\ ■ /nn\ 

= kEj(if 2 -if 2 m, (22) 

where we have used that d V2 v(ipi, f 2 ) vanishes at the 
minima t^™ 1 "'™ and only the term originating from the 
coupling term in v(ip±, f2) is relevant. Since ipi is at most 
of order u>o,i, we obtain the estimate dt£%o ~ ^kujQ^Ej. 
The smallness of plz, guaranteeing adiabatic motion of 
the classical junction (this requires a correspondingly 
large C\) then follows from 

The relaxation of the quantum degree of freedom f2 is 
another important element in our discussion. The dis- 
sipation described by the normal resistance R in the 



RCSJ-equation of motion, cf. Eq. (fTJ), leads to typical 
finite lifetimes r 2 ~ RC2 = l/2a2W p ,2 of the excited 
states \l) of the quantum junction^ A much longer life- 
time T2 shows up if the quantum junction is trapped in 
a local ground state of a side- well in v e s(f2), cf. Fig. 
[5] The decay then is protected through a large barrier 
Et>, enhancing the typical lifetime to a value T2, with 
r 2 <C f 2 oc exp(2 r yy / Ej /AE C ,2)^ Thus, there are two 
(extreme) ways how a highly exited state in the 'super- 
well' of v e g(if2) can decay, see Fig. [5J either via states 
within the superwell involving the typical lifetime r 2 , or 
via states involving a tunneling process, resulting in an 
exponentially larger decay time of the order of T2 . 
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FIG. 8: Illustration of decay sequences of an excited 
'superwell'-state of the quantum junction in the effective po- 
tential Vcs (952), Eq. Q. Solid arrows indicate a decay se- 
quence involving only states within the 'superwell' with a 
typical decay time t%. The dashed arrows illustrate another 
(extreme) decay sequence, involving localized states in the 
side-minima of v c g{ifi2)- The decay time then involves the 
typical lifetime 2> T2 of local ground states within side- 
wells. 



In the following, we investigate three regimes and de- 
rive the corresponding effective critical current j c (k). We 
begin with the case of adiabatic ramping (Sec. IIV A[) . 
where we assume that the ramping is slower than the 
typical relaxation time of the heavy junction and the 
timescales between (p% and if 2 separate; i.e., inequality 
(|2"3"|) holds. We show that the critical current j c (k) ob- 
tained in the dissipative case, Sec. IIII1 is only slightly 
altered due to the (small) finite amount of kinetic energy 
the classical junction can store. Next, we discuss the case 
of fast ramping fSec. HVFlj) . where the phase <f \ remains 
effectively 'frozen' during the current ramping and quite 
an appreciable amount of potential energy is converted 
to kinetic energy in the motion of the classical junction, 
leading to a pronounced reduction of j c (k). Both discus- 
sions require a very large value of Ci/C 2 —> 00 to guar- 
antee the absence of Landau-Zener tunneling due to the 
motion of ipi. Subsequently, we discuss the consequences 
of a moderate ratio C\/C 2 < fO 4 accessible in today's 
experiments in the last part of this section (Sec. IIV C[) . 
Then, if 2 can remain trapped in a localized state in a side 
well during the motion and liberation of the heavy junc- 
tion. The probabilistic effects emerging for this situation 
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change the nature of the phase diagram qualitatively. 

A. Adiabatic ramping of the bias current 

We start from the system with fixed coupling k at j = 
and consider a state localized at (pi = cp2 = 0. We assume 
slow current ramping on the dissipative timescale of the 
classical junction's motion, 

dtj « 1/n, (24) 

where t\ ~ RC\ denotes the classical relaxation time of 
the heavy junction. The inequalities (|24[) (slow ramping) 
and (|23p (slow motion of ^i) guarantee that the quantum 
junction remains in the global ground state, while the 
classical junction follows its local ground state near ip\ = 
arcsin j . 

The effective potential for the heavy junction then is 
determined by the ground state energy So(ipi). In our 
estimate 

e Q { Vl )^Ejv[ Vl ^t h { Vl )] (25) 

we neglect the correction due to the ground state energy 
Eq hujp^/'i (as measured from the bottom of the po- 
tential) of the light junction^ the phase ip^ oh (ipi) refers 
to the global minimum of the light junction. Further- 
more, we ignore the small splitting 2A2 at the avoided 
level crossing (cf. Fig. [71(a)), replacing it with a sharp 
kink. 

With the above approximations, the force acting on 
the heavy junction originates from the classical force ex- 
erted by the quantum junction residing in the global mi- 
nimum. The situation is thus the same as in the strong 
damping regime: the 'flux-entry lines' j^ n (k), marking 
the current where the light junction is allowed to tunnel, 
are still given by Eq. (fT2| . When the classical junction 
has relaxed to the bottom of a minimum, that minimum 
has to disappear in order for ipi to enter the running 
state. Given the expression (125|) for eo(ipi), the local 
minimum in £o(</?i) corresponding to the n-th side mi- 
nimum in v c ff (<^i, ^2) turns into an inflection point at 
j = jc n (k), Eq. (fT4| . This condition again agrees with 
the one for strong damping. 

A difference to the previous dissipative situation is 
given by the possibility to transform potential into ki- 
netic energy. Ramping j at fixed k past the flux-entry 
line jc n {k) triggers a phase slip in ip 2 and deforms the 
potential well in eo( < / 9 i) for the classical junction. The 
latter finds itself on the slope of the well and gains ki- 
netic energy while sliding down towards the new mini- 
mum. The classical junction then may enter a running 
state if this kinetic energy gain is sufficient to surpass 
the barrier blocking the newly formed (n-th) minimum. 
Different from the overdamped case, the system then ex- 
hibits hysteretic behavior, as illustrated in the dynamical 
phase diagram Fig. [HI where the running state can be en- 
tered along extensions of the j^ n (k) line away from the 
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FIG. 9: Phase diagram of the dynamically asymmetric dc- 
SQUID for weak damping and adiabatic ramping, (dtj <C 
Tj -1 ). The critical line j c (k) is nearly identical with the one in 
the strong damping case, but is hysteretic along extensions of 
the critical lines j£ n (k). These are due to the nonzero kinetic 
energy acquired by the classical junction when sliding down in 
the new well formed after the tunneling of <fi2 at j^„(k) and 
allowing tpi to surmount a residual potential barrier in the 
effective potential eo(fi)- The inset illustrates the condition 
determining the termination point of the hysteretic line along 
j^x(k). The three solid dots mark points of equal potential 
energies; the energy gained by the classical junction due to the 
phase slip of ip2 just suffices to climb the top of the remaining 
barrier along ipi. 



phase boundary. The point where these extensions ter- 
minate is determined by the condition that the energy 
right after the tunneling of ip 2 is equal to the energy at 
the top of the barrier, as illustrated in the inset of Fig. 
[HI If the coupling k is too small to lower the barrier suf- 
ficiently (cf. Fig. [7J (b)), the state of the heavy junction 
upon reaching j = j^ n (k) is transformed into a localized 
excited state of the newly formed potential well in Eo(ipi). 
Given the slow ramping d t j -C rf , the heavy junction 
relaxes to the ground state of the well and a further in- 
crease in j is necessary to remove the barrier completely, 
hence the critical line jumps forward to j c (k) = 3cn(k)- 

B. Fast ramping of the bias current 

We now turn to the case dtj 3> w P) i, where the ramp- 
ing of j is instantaneous with respect to the (massive) 
dynamics of the heavy junction; the final current j has 
to remain below ~ 0.7, since fast ramping beyond this 
value allows the classical junction to overcome the po- 
tential barriers through conversion of potential to ki- 
netic energy. Furthermore, we study the setting where 
cjp i <C T2 ; the quantum junction then has relaxed 

to the ground state (at ip 2 nm ' n ) in the effective potential 
feff (^2) before any motion of the heavy junction sets in. 
Hence, initially, the effective potential for the classical 
junction is given by the ground state energy eo(<pi). The 
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absence of Landau-Zener tunneling (due to the inequality 
([2"3"| ) then assures that the effective potential is given by 
So((pi) during the entire motion of the classical junction. 

We then can find the criterion for <*p\ to enter a run- 
ning state: During the fast ramping of the current j, the 
heavy junction remains frozen at tpi = 0. We have to 
check whether the potential energy of the heavy junction 
is sufficient to overcome all barriers in £o(</?i > 0); i.e., 
we have to inspect if the maximum of So((pi) in the in- 
terval < ipi < 2ir is realized at ipi = 0. As soon as this 
condition applies, the classical junction becomes delocal- 
ized. This analysis has to be performed numerically, and 
the result is displayed in Fig. [lOl As a first overall re- 
sult, we note that the critical current is shifted to lower 
values due to the instantaneous change of the effective 
potential for tpi by the fast ramping of the current and 
by the quantum decay of (f2, allowing the heavy junction 
to transform potential into kinetic energy. 

Another prominent change is in the form of the tran- 
sition line, specifically, the 'cut' of the first tip of j c {k): 
The first segment at large values of k, cf. Fig. [lOl cor- 
responds to the conventional case for a + -type line; the 
criterion for derealization involves the full conversion 
of potential to kinetic energy of the heavy junction and 
the subsequent tunneling of the quantum junction right 
at the turning point of the classical junction, cf. inset 
(1) of Fig. [TUJ Following the j+(fc)-line further with de- 
creasing j, we reach a point where the potential energy 
of the classical junction after the phase slip event is no 
longer sufficient to overcome the next barrier and an ad- 
ditional increase in j is required to provide the missing 
energy. The segment (2) then involves partial conversion 
of the initial potential energy within the two wells before 
and after the phase slip event. Finally, moving along the 
segment (2) towards lower k the value of ipi where the 
phase slip occurs decreases and finally reaches tpi = 
(cf. Fig. [THl inset (3)). The critical line then continues 
along the third segment which is of the typical ~-type, 
with the quantum junction undergoing tunneling before 
the classical junction starts moving. After tunneling, a 
large current is required to lower the barrier in the cosine 
potential of the classical junction, cf. Fig. [TQl inset (3). 
Note that the differentiation between positively sloped 
segments of type (1) and type (2) persists to segments 
determined by a higher fluxon index n, although the dif- 
ference is invisible on the scale of Fig. [TU1 



C. Moderate ratios of C1/C2: Probabilistic effects 

Above, we have assumed the limit of a very large ratio 
C1/C2. In the experiment, the dynamical asymmetry of 
the SQUID is implemented by shunting one of the two 
junctions with an external capacitance. 33 This technique 
allows to reach large capacitance ratios of the order of 
C1/C2 S 10 4 - Given this finite mass ratio, one may en- 
ter a new dynamical regime where the classical junction is 
pulled out of its metastable state more efficiently. This is 
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FIG. 10: Effective critical current j c (k) (solid) of the current- 
biased dc-SQUID with strongly asymmetric parameters as a 
function of bias current j and inverse inductance k oc 1/L. 
We assume instantaneous ramping of the bias j and a strictly 
classical and adiabatically slow (as compared to the dynam- 
ics of the quantum junction) motion of the weakly damped 
heavy degree of freedom ipi. For comparison, the critical line 
of the dissipative case is also displayed (dashed line). The 
inset illustrates the three segments constituting the first part 
of j c (k), where solid dots indicate points of equal potential 
energy. Inset (3) shows the situation right at the tip between 
segments (2) and (3). 



realized when the quantum junction attains a high energy 
state residing at the opposite side of the parabolic poten- 
tial (large values of (f2 in Fig. [5]) and holds on to it while 
dragging the heavy junction out of its metastable state. 
This requires a fast^l current ramping d t j 3> I/T2 (on 
the scale of the intervalley motion of the light junction), 
pushing the quantum junction to high potential energies 
during the ramp-up process. Furthermore, the condition 
w p,i ^ 1/^2 has to be satisfied to allow for motion of the 
heavy junction before relaxation of the quantum junction 
(from a side-well local ground state) back to the global 
minimum. 

In order to simplify the situation, we assume a (weak) 
finite damping of the classical junction, u)~\ <C T\ < T2, 
assuring that the classical junction relaxes before tunnel- 
ing of the quantum junction (this assumption is merely 
a convenience, allowing us to ignore the availability of 
kinetic energy for the classical junction upon fast ramp- 
ing). We then start from the local ground state of the 
original potential well at ipi = (p% = arcsinj. The decay 
proceeds in the direction of ip2 in the effective potential 
v c ff [arcsinj] ((^2) and can end up in any of the lower local 
ground states for the quantum junction. We denote the 
energy of the local ground state of the quantum junc- 
tion in the n-th well by £q (not to be confused with e n , 
the n-th excited state in the spectrum of the quantum 
junction). We estimate these energies as 



~n { \ { mm,n\ 

e (^i) w v(tpi.,ip 2 ' ), 



(26) 
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where <£>™ ln ' n (c/?i) is the solution to Eq. (|8|) near if2 ~ 2im 
at fixed tpi and we have neglected the zero point energy 

The dynamics of ipi depends strongly on the fluxon 
index n of the state of the quantum junction after the 
decay: the larger n, the stronger is the force on ipi fa- 
cilitating its escape. We can define two extreme values 
for the critical current, a lower limit (k) arising from 
the largest attainable fluxon index n, and an upper limit 
j {2) {k) associated with the index n of the global mini- 
mum of the quantum junction. We estimate the maxi- 
mal index n by comparing the energy of the state of the 
quantum junction before its inital decay with the height 
of the barrier blocking the access to the n-th minimum; 
the numerical solution of the relation 

u e ff[yi = arcsinj]((^ 2 = arcsinj) = v cS (ip™ ax ' n ), (27) 

with ( y 3 I 2 nax '™ the solution to Eq. © near ipf ax w (2n - 
l)ir at if i = arcsinj provides us with the line j^^~(k). 
The subsequent stability analysis of <pi generates the lines 
jc~n~ V s ) i their combination into j^ 1 ' (k) is displayed in Fig. 

ni 
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FIG. 11: Effective critical currents of the dynamically asym- 
metric SQUID for T2 S> tdpi and fast ramping. The system 
always turns resistive for j > j^ 2) (k) but never turns resistive 
for j < jc(k). For j^ 1 ' < j < ji 2 ' 1 the derealization of the 
system is determined by the statistical nature of the decay 
process of the quantum junction involving the side minima. 



The shape of j {2) (k) is found as described before: We 
determine the index n of the global minimum of w e ff (^2) 
and deduce the motion of ipi in the effective potential 
£q (tpi). The reduction in critical current due to the mas- 
sive dynamics of <pi is not as pronounced as in Fig. QJ3 
since the maximal critical current j ( c 2) (k) involves the pro- 
cess where the classical junction dissipates most of its ki- 
netic energy: the quantum junction decays to the (n — 1) 
fluxon state, the classical junction relaxes to its local mi- 
nimum near ipi > arcsinj, whereupon the quantum junc- 
tion decays further to its global ground state with fluxon 
index n. Only the potential energy gained by the classi- 
cal junction in the last step can be transformed to kinetic 
energy, making the resulting line shown in Fig. II II differ 
from j c (k) in the overdamped case (Fig. [5]). 



The two critical lines jc } (k) and j {2) {k) define a 
broad intermediate regime, see the grey area in Fig. 
[Til where the decay of the system is of probabilistic 
nature: whereas no Munchhausen decay is possible for 
j < jc 1} (k)i the decay of the system in the intermediate 
regime depends on which minimum the quantum junc- 
tion cf2 decays to. For j > j {2) {k) the system is always 
bound to decay. 

The critical line j {2) (k) crosses twice the critical line 
j^(k) near k ~ 0.5 and for n = 2. This special situation 
arises when the most distant accessible minimum and the 
global minimum are neighbors (with indices n = 2 and 
n = 1) and hence strongly interrelated. With decreas- 
ing k, both the height of the maximum near tp 2 ~ 2>tt 
(determining j^ik)) and the position of the minimum 
at if2 ~ 47r (n — 2) with respect to the one at if2 ~ 2?t 
(ri = 1; the competition between these minima defines 
j {2) {k)) are lowered. The latter effect is slightly stronger, 
leading to the crossing of the lines under a small angle. 
Upon further lowering of k, the relevant fluxon index for 
jc } (k) increases twice as fast as for j {2) (k), preventing 
any further crossings of the critcal lines. 

To guarantee reasonable time scales in the experiment 
requires a fast quantum dynamics near the global mini- 
mum, hence the ratio Ej/fku Pi 2 should be small, of order 
unity. For a fast ramping, this is in conflict with the re- 
quirement of large T2 in all side- wells; as a result, not all 
reachable side minima may be relevant for the dynamics 
of ipi (i.e., T2 oj~\ does not hold), effectively shifting 
the lower critical current j^, 1 ' (k) to larger values. 

So far we have neglected any quantum corrections for 
the variable ip\\ however, with a finite ratio C1/C2 such 
corrections may become relevant. Modifications due to 
quantum fluctuations in the heavy degree of freedom are 
twofold: On the one hand, we have disregarded the tun- 
neling of ifi through a small residual barrier, leading to 
a broadening of all lines of the "-type. On the other 
hand, we have assumed the dynamics of ip 2 to occur in 
the effective potential v e s[tpi](ip2) , with tpi fixed, while 
for a finite value of C\ the ground state wave function of 
the heavy junction acquires a finite width of the order of 
Sipi ~ (SEd/Ej) 1 / 4 . This leads to a broadening of all 
lines of the + -type, which can be interpreted as tpi taking 
part in the tunneling of <f2- 



V. EXPERIMENTAL IMPLEMENTATION 

The running state following a Munchhausen decay is 
associated with a finite time-averaged voltage across the 
device. The above results thus can be tested experimen- 
tally by measuring the voltage drop as a function of the 
applied bias current J and of the inductance L of the 
SQUID. 

Mapping any part of j c (k) requires changing the cou- 
pling k oc L~ l . The latter is determined by the total 
inductance of the loop L = L geo + Lkm, comprising the 
geometric (L goo ) and kinetic (Lkin) inductances of the 



11 



loop. The former is of the order of the loop perime- 
ter and accounts for the magnetic field energy. The ki- 
netic inductance relates to the geometric quantity via 
^kin ~ L gco X 2 jr 2 , with r the radius of the wire and A 
denoting the superconducting penetration depth. Usu- 
ally, the inductance is dominated by its geometric part. 
However, a large inductance, as desired for our setup, 
can be installed by using a thin wire. Furthermore, the 
inductance is fixed after fabrication and its modification 
(e.g., via installing a diamagnetic shield) is rather dif- 
ficult. Another way to trace the critical current line is 
obtained by applying an external magnetic flux $ e to 
the sample, in which case the potential Eq. ([2]) has to be 
replaced by 

v((fii, if 2) = 1 - cost^i + 1 — cos<p 2 
k 

~ j(<P! + <p 2 ) + -((fi-<P2- 27T$ e /$ ) 2 - (28) 

The critical current j c then can be studied as a function 
of $ e - The analysis proceeds in the same way as before: a 
change in k altering the opening angle of the parabola in 
v((fi, LP2) is replaced by a shift in the parabola's position 
due to the applied flux <I> e . 
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FIG. 12: Phase diagram of the dynamically asymmetric 
SQUID in the strong damping case as a function of the exter- 
nally applied flux $ e and bias current j. Here, the coupling 
k cc L~ x is fixed at a value k = 0.06; the critical lines in the 
diagram correspond to n — 2. The critical line is periodic in 
<& e with a period $0, a consequence of the invariance under 
the replacements $ e /$o — > ^e/^o ± 1 and n — > (n ± 1). 



As an illustration, we discuss here the most simple case 
of an overdamped dynamics, analogous to Sec. IIII1 In a 
first step, the system relaxes to the ground state at the 
minimum of the potential v((pi,ip2), Eq. ([28]) , as given 
by the solutions <p™ n , </2™ ln of 

sin^f") -j + k (Vf 1 ~ Vf n - 2tt|^ = 0, (29) 

sin(^ in ) -J-k (pf n - <fif n - 2tt|^ = 0, (30) 

near cpf in w ip 2 nin ks 0. The critical line j+„(<£> e ) then 
derives from the condition 



where both the shift of the parabola v e g (^2) by 27r$ e / $0 
and the flux dependence of ip™ ln are accounted for. Sim- 
ilarly, Eq. ([14]) for k~ n (j) is changed to 
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(2n - l/2)?r+ arcsin(2j - 1) + 27r$ e /$ ' 



(32) 
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(2n - 1)tt - (^r n (*e) + 27T$ e /$ ' 



(31) 



from which we find the critical lines j~ n (Q e ). The results 
of Eqs. ([29]) - ([32]) are invariant under the shifts $ e /®o — » 
Q e /$o ± 1 and n — > (n ± 1). The Miinchhausen decay 
thus can be studied within a finite interval, e.g. $> e /$o G 
[—0.5,0.5], and the resulting phase diagram is displayed 
in Fig. [H 

Biasing the SQUID with an external magnetic flux has 
a convenient side effect: For negative flux $ e , the crossing 
point of jcn{k) an d 3cn(k) i s shifted to larger values of 
k and the tunneling barrier for ip 2 is lowered, allowing 
to study the Munchhausen effect using a SQUID with 
considerably smaller inductance and hence smaller size. 



VI. CONCLUSION 

We have shown how a system consisting of a 'quan- 
tum' degree of freedom coupled to a 'classical' one, im- 
plemented with a two-junction SQUID, can perform a 
complex decay process out of a zero- voltage state involv- 
ing quantum tunneling of the 'light' junction and classi- 
cal motion of the 'heavy' junction. We have analyzed the 
cases of strong and weak damping, for different methods 
of preparation (fast and slow ramping), and for various 
ratios of the capacitances C\ / C2 and have found the ef- 
fective critical current j c (k) exhibiting a common char- 
acteristic shape. Its zig-zag like structure originates in 
two competing effects: while the number n of flux units 
which can enter the SQUID is increased with decreasing 
k, the current redirected over the heavy junction at given 
n is decreased. For junctions with equal critical currents, 
the critical line approaches the value j c (k — > 0) — » 1/2 
in the dissipative case; asymmetries in the critical cur- 
rents of the two junctions lead to a shift of the critical 
line j c (k) to values straddling the critical current of the 
heavy junction. Going over from dissipative to massive 
dynamics, the kinetic energy stored in the system allows 
to better overcome the potential barriers and the critical 
line shifts to smaller values of j; furthermore, hysteretic 
effects, similar to the situation encountered for under- 
damped junctions, and even probabilistic behavior may 
show up. 
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